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Stimulated by the small/large Fermi surface controversy in the cuprates we consider a small num- 
ber of holes injected into the bilayer antiferromagnet. The system has an 0(3) quantum critical point 
(QCP) separating the magnetically ordered and the magnetically disordered phases. We demon- 
strate that nearly critical quantum magnetic fluctuations can change the Fermi surface topology 
and also lead to spin charge separation (SCS) in two dimensions. We demonstrate that in the phys- 
ically interesting regime there is a magnetically driven Lifshitz point (LP) inside the magnetically 
disordered phase. At the LP the topology of the hole Fermi surface is changed. The position of the 
LP, while being close to the position of the QCP is generally different. Dependent on the additional 
hole hopping integrals t' and t", the LP can be located either in the magnetically ordered phase 
and/or in the magnetically disordered phase. We also demonstrate that in this regime the hole spin 
and charge necessarily separate when approaching the QCP. The considered model sheds light on 
generic problems concerning the physics of the cuprates. 

PACS numbers: 74.40.Kb, 74.72.Gh, 75.10.Jm, 75.50.Ee 



I. INTRODUCTION 

It is well established that electron/hole diffraction from 
a static magnetic order in a conductor can influence the 
FS topology; this is FS reconstruction. In addition it 
is also well established that spin and charge are sepa- 
rate in one-dimensional systems^. This results in two 
interesting but important generic questions; (i) Can dy- 
namic magnetic fluctuations also drive a change in the 
FS topology? (ii) Is spin-charge separation possible in 
two-dimensions, and if so what is the meaning of the sep- 
aration. In the present work we address these two generic 
problems and demonstrate that these two problems are 
remarkably related. These two important issues have re- 
cently attracted alot of attention due to their close tie to 
the high T c superconductivity in the cuprates. 

To address the two generic problems of interest we con- 
sider a small number of holes injected into the bilayer 
Heisenberg antiferromagnet where the interlayer coupling 
drives the strength of the magnetic fluctuations. The 
system has an 0(3) quantum critical point (QCP) sep- 
arating the magnetically ordered and the magnetically 
disordered phases. We show that indeed purely dynamic 
short range AF correlations in the absence of a static AF 
order can cause a LP in which the topology of the FS 
is changed. The position of the LP, while being close to 
the position of the QCP is generally different. Depen- 
dent on the additional hole hopping integrals t' and t" ', 
the LP can be located either in the magnetically ordered 
and/or in the magnetically disordered phase. Vojta and 
Becker 3 also observed a LP in a similar model. However, 
in their study the LP was always in the AF ordered phase 
and therefore the central issue of the dynamic magnetic 
fluctuation driven LP has not yet been addressed. 

We also demonstrate that the hole spin and charge nec- 



essarily separate when approaching the magnetic 0(3) 
QCP when the LP is in the magnetically disordered 
phase. The possibility of spin-charge separation (SCS) 
in 2D has been previously studied using the slave-boson 
method^—. However, the slave-boson method as applied 
to the t — J model implies the SCS ad hooi. This is not 
the case in the present work, therefore, the precise mean- 
ing of SCS in the present work is different from that of 
the slave-boson method. 

Our interest in the two generic problems is motivated 
by the cuprates, which, in our opinion, manifest both the 
LP as well as SCS. It is therefore appropriate to explain 
the connection with the cuprates which we do below. A 
reader who is interested in these generic problems but 
not necessarily their application to the cuprates can go 
directly to Eq. (1) where we start the analysis. 

One of the central issues in understanding high-T c su- 
perconductivity is whether it originates from a Fermi 
liquid or from a Mott insulator. While there is no 
consensus on the problem, there are experimental in- 
dications including angle-resolved photoemission spec- 
troscopy (ARPES) that the transition from small to large 
Fermi surface occurs in the hole doping range 0.1 < x < 
0.1 Fermi arcs^ observed below this doping level 
have triggered various theoretical studies and models. 
However, recent ARPES data supports the existence of 
a shadow band that completes the arc into a pocke t 13 ! 14 . 
The existence of hole pockets is in agreement with several 
calculations which considered dilute holes dressed by spin 
fluctuations, based on doping a Mott insulator——. On 
the other hand, a large Fermi surface as expected from a 
Fermi liquid approach is observed in ARPES studies in 
the optimally to over doped cuprates. Naturally, this im- 
plies that there is at least one topological Lifshitz point 22 
(LP) in the doping range 0.1 < x < 0.15, where the Fermi 
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surface changes from small to large in contradiction to 
the Luttinger theorem. A phenomenological description 
based on Fermi liquid picture^ 3 -, as well as dynamical 
mean-field theory calculations within Hubbard model 2 ^ 
have been proposed to describe this LP. 

Magnetic quantum oscillation (MQO) in underdoped 
YBa2Cu30^ supports the small pocket scenario, in con- 
trast to the large FS observed on the overdoped side 2 ^. 
The above observations also suggests the existence of 
a LP at which the topology of the FS changes. The 
MQO measurements were performed in very strong mag- 
netic fields, up to 80T. It has been suggested 2 ^ that the 
field induces a static magnetic structure and the struc- 
ture causes the small Fermi surface reconstruction. On 
the one hand, the MQO experiments suggest the small 
FS was observed up to 12% doping 2 ^, and it is unlikely 
that even an 80T field can generate a static AF order 
at such high doping. On the other hand the short range 
dynamic AF correlations always exist in the cuprates; 
this has been supported by recent RIXS measurements^ 
which have remarkably demonstrated that such correla- 
tions are practically doping independent, from Mott in- 
sulator to optimal doping. Based on this data one can 
conjecture that the cuprates are always close to magnetic 
criticality. This motivates us to study if the LP can be 
driven by short range, purely dynamic AF correlations. 
We consider a bilayer model for the sake of performing 
a controlled calculation. However, we believe that con- 
ceptually our conclusions are equally applicable to both 
single and multi-layer cuprates. 

Besides the possible topological transition, this partic- 
ular doping range 0.1 < x < 0.15 has also been associated 
with a quantum critical point (QCP) where the static 
magnetic ordering becomes fully dynamic. Below this 
doping range, neutron scattering indicates a commen- 
surate or incommensurate magnetic ordering depending 
on the doping level as well as the coupling between the 
Cu02 planes 3 - -"— . The static magnetic order vanishes 
within this doping range, as also confirmed by nuclear 
magnetic resonance (NMR) wipeout j 34 i 35 and recently 
confirmed by muon spin rotation (/iSR) data 3 -^. The re- 
duction of the static magnetic moment has been calcu- 
lated within a low energy effective theory 3 ^, and a QCP 
around x ~ 0.11 is predicted, close to the experimen- 
tal value x w 0.09 observed in YBa 2 Cu 3 ?/ ^ 3 -^^. It 
is then intriguing to ask whether the topological LP is 
related to this magnetic QCP, especially how the change 
of magnetic excitations can influence the electronic prop- 
erties. However, currently there is no precise measure to 
discern whether the topological transition takes place in 
the magnetically ordered or the fully dynamic phase. 

Now we turn to the discussion of SCS and its relation 
to magnetic quantum criticality. In La2-^Sr :c Cu04 the 
QCP is smeared out because of disorder. However, in 
YBa2Cu30 2y the QCP is located experimentally at dop- 
ing x « 0.09 (y « 6.47)^^^. SCS manifests itself 
in the cuprates as a phenomenon related to the mag- 
netic QCP mentioned above. In the magnetically ordered 



phase, the pseudospin that marks the sublattice is a good 
quantum number for holes. The meaning of partial SCS 
is then the unusual way the pseudospin interacts with 
an external magnetic fiel d 40 ! 41 and this is the meaning 
of the partial SCS in the magnetically ordered phased. 
On the other side of the QCP, the lack of long range 
order indicates that spin and charge are recombined, as 
expected from the Fermi liquid approach. However, it is 
unclear how SCS develops as one approaches the QCP 
from the magnetically disordered phase, or stating from 
another aspect, if there is a certain residual effect from 
SCS that affects the hole motion in the magnetically dis- 
ordered phase. In addition, since the QCP and LP are 
located in close vincinity, is SCS related to or participate 
in the mechanism that causes FS reconstruction. In the 
present work we analyse the process of SCS at the QCP. 
The model considered here has only commensurate mag- 
netic ordering, so we put aside incommensurability in the 
cuprates. 

Our recent work 42 suggests that the two seemingly un- 
related issues of FS reconstruction and SCS may origi- 
nate from the same mechanism, namely nearly critical 
quantum fluctuations. By studying the bilayer AFM 
in the magnetically disordered phase, we show that dy- 
namic magnetic fluctuations, driven by the ratio of the 
interlayer to the intralayer coupling J_l/J, has a dra- 
matic effect on both the dispersion and spectral function 
of a single hole. The system under consideration has a 
0(3) QCP at (J_l/J) « 2.31, close to which the disper- 
sion has minima at k = (±7r/2, ±7r/2). This results in 
small hole pockets similar to that in the cuprates at small 
doping. The pockets are formed due to strong in-plane 
AF correlations which diminish the nearest site hopping 
t. Upon increasing J± the in-plane AF correlations are 
reduced and the dispersion minima gradually shift to- 
ward k = (±7r, ±7r) reaching this position at some value 
J± > Jf p . This is the position of the LP where the 
shape of the single hole dispersion changes. In-plane AF 
fluctuations are negligible at very large J_l, so the hole 
dispersion recovers the shape expected from a Fermi liq- 
uid approach, but the bandwidth is a factor of two re- 
duced due to J±. Note that we examine only the single 
hole case to demonstrate the mechanism, so the notion 
of a FS is ambiguious. Nevertheless, it is known that the 
SCBA typically converges at small doping up to x ~ 0.1 
or so, and the rigid band approximation in most mod- 
els is valid. Therefore we expect our results to be valid 
at small doping, where the FS changes from four small 
pockets centered around k = (7r/2,7r/2) at J±_ < Jf p to 
one small pocket centered at k = (jr, it) at J± > J^ p '. 
Although it is unclear at present if this model can be 
realized in certain materials and if the driving parame- 
ter J±/J can be controlled externally, one clear predic- 
tion is that if an MQO experiment is available for such a 
compound, then one would observe the MQO frequency 
changes by a factor of four across the LP. 

Taking these issues back to the cuprates, an intrigu- 
ing argument for the relevance of our calculations comes 
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from recent resonant inelastic x-ray scattering (RIXS) 
experiments^. RIXS indicates that magnetic fluctua- 
tions exist even in the overdoped regime, and is prac- 
tically unrenormalized as in the undoped parent com- 
pound. Therefore purely dynamic magnetic fluctuations 
are very likely to be the driving force behind FS recon- 
struction, even though the LP in the cuprates takes place 
in a doping region that has no static magnetic order. 
However, one has to keep in mind that the magnetic 
fluctuations in the cuprates, regardless of the possible 
multilayer structure in a unit cell, are likely to be purely 
2D, and the interlayer coupling is unlikely to be a crucial 
driving force for the magnetic fluctuations. Therefore, 
despite the appealing conclusions from the bilayer AFM, 
a theory that can address the 2D magnetic fluctuations 
in the cuprates will be much more challenging, and is yet 
to be formulated. 

The purpose of this extended version is to provide a 
detailed formalism for the results presented in Ref. E2I 
and to clarify several points that we omitted to address. 
We provide details for the bond operator formalism that 
describes the magnetic fluctuations in terms of triplon ex- 
citations. This mean field formalism yields a good agree- 
ment with Quantum Monte Carlo and series expansion 
result s 46 ! 47 ! 56 , and has the advantage that hole dynamics 
can be treated within a simple self-consistent Born ap- 
proximation(SCBA). As mentioned in Ref. 42, the hole- 
triplon vertex has a discrepency with a previous study 43 , 
so we provide detailed calculations to justify our vertex. 
To support the argument of magnetic fluctuation induced 
FS reconstruction, we have drawn evidence from series 
expansion calculations to justify our SCBA approach. In 
this long version we provide numerical results of the se- 
ries expansion calculations and make a direct comparison 
with the SCBA calculations. 

The structure of the paper is the following. In Sec. [TT] 
we formulate the model and discuss the relevant set of 
parameters. In Sec. IIII1 we present a brief overview of 
the bond operator formalism for the double layer Heisen- 
berg model and reproduce the results of Refs . 44 i 45 for 
spin excitations in the the magnetically disordered phase. 
In Sec. |TVl we derive the hole- triplon verticies for the 
quantum disordered phase, and compare them with pre- 
viously known expressions. In Sec. [Vlwe introduce the 
self-consistent Born approximation (SCBA) and calcu- 
late the hole Green's and spectral functions. Details of 
the numerical calculations are also presented. In Sec. [VI] 
we outline the series expansion methods used in our study 
and in Sec. IVIII discuss the results of our numerical cal- 
culations and comparison with dimer series expansions. 
In Sec. IVIIII we address the issue of SCS in the vicinity of 
the QCP. Finally, in Sec. HXI we present our conclusions. 



II. MODEL, PARAMETERS 

We consider the t — t f — t" — J — J± model on a bilayer 
square lattice at zero temperature 




FIG. 1: (color online) Schematic diagram of the Heisenberg 
bilayer antiferromagnet doped with a single hole. We show the 
hole hopping integral parameters Uj, the antiferromagnetic 
exchange J between spins on neighbouring sites as well as the 
interlayer antiferromagnetic exchange J±. 



H = J^(S«.sf) + sf).sf) + J^sW.sf) 

(ij) i 

~ ^1 ( c ila C jlo- + Cj 1(J Cn a ) 

(hj) 

= Hj^ J± + H t ,t f ,t" (1) 

where c\ a is the creation operator of an electron with spin 
a =t, I at site i on the top plane, = \c\ cf '^ v Ci V , and 
Uj = {t,t',t"} is the hopping integral between nearest, 
next-nearest, and next-next nearest neighbour sites re- 
spectively as shown in Fig. [TJ The superscripts (1), (2) 
in (pQ) indicate the layers. To compare parameters on a 
similar energy scale hereafter we set J = 1. A no-double- 
occupancy constraint is imposed on the system. This 
means that at half filling (one electron per site) there 
are no mobile charge carriers and we have a Mott insu- 
lator. The Hj j j ± term describes the antiferromagnetic 
coupling in each layer as well as between the two layers. 
It is well known that without holes (half-filling) the model 
has an 0(3) magnetic QCP at J± = 2.525 46 i 47 separating 
the magnetically disordered and the AF ordered phases. 
The hopping integrals t,t' ,t" result in charge dynamics 
if a hole is injected into the system. Note that hopping is 
only allowed in the top plane, so a mobile hole can only 
be injected into the top layer. The longer range hopping 
integrals t' ,t" are crucial as we will explain later. Note 
also that since we consider the zero temperature case, the 
magnetic ordering in the AF phase is consistent with the 
Mermin- Wagner theorem. 

A similar extended t — t' — t" — J model has been widely 
applied to study the magnetic properties of the prototype 
bilayer cuprate YBa2Cu306+ 2/ in the regime J±/J <C 1, 
as indicated by fitting the magnon dispersion via neutron 
scattering dat a 48 i 49 . In the current study, we investigate 
the magnetically disordered phase in the opposite limit, 
J±/J> 1, as addressed below. 

In the present study we focus on small doping, 
such that the magnetic fluctuations are not influenced 
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t/J 


# / / T 
t 1 J 


iff I T 
t 1 J 


3.1 


-0.8 


0.7 


3.1 


0.0 


0.0 


0.5 


0.0 


0.0 


0.5 


0.0 


0.1 



TABLE I: Hopping integral parameters tij in units of the an- 
tiferromagnetic exchange J. Here we consider both "strong" 
(t/J = 3.1) and "weak" (t/J = 0.5) coupling. 



by doping. Therefore we can assume that the magnetic 
fluctuations and the QCP is driven by the interlayer cou- 
pling J± alone. The holes fill the rigid band formed by 
the magnetic quantum fluctuations. To address the prob- 
lems formulated above it is sufficient to calculate the sin- 
gle hole Green's function. Certainly for sufficiently high 
doping the holes start to influence the magnetic fluctua- 
tions and hence the rigid band approach fails. However, 
it is not necessary to go to such high doping to draw our 
conclusions. Such an approach is only possible because 
the magnetic dynamics are driven by J± and are indepen- 
dent of the hole concentration. This is a significant sim- 
plification compared to the t — J/Hubbard model, where 
doping is the only "handle" . 

Although the current study does not intend to simulate 
the cuprates directly, we choose the parameters relevant 
to those in the cuprates such that the interplay of various 
mechanisms can be compared on a similar energy scale. 
The in-plane antiferromagnetic coupling in the cuprates 
as measured by two magnon Raman scatterin g 50 ! 51 is 
J = 125 meV. As will be addressed in later sections, 
we investigate several different hopping parameters and 
show the importance of t' and t" . To make a relevant 
connection to the cuprates, we choose one of the pa- 
rameters close to that of YBa2Cu307, which in the first 
principle study yields t = 386 meV, t' = —105 meV, 
t" = 86 meV^. In order to compare the interplay of var- 
ious mechanisms on a similar energy scale we set J = 1. 
These values have also been applied to the analysis of 
photoemission in S^CuC^C^r^ an d yield a good agree- 
ment with the ARPES data. In the current approach we 
use the values provided in Table. [U 

The motivation of using these parameters is that they 
cover both "strong" and "weak" coupling regimes. Here 
we denote "strong" with large t and "weak" with small 
t. The importance of the "strong" coupling regime is 
stressed by its relevance to the cuprates. The "weak" 
coupling regime is particularly useful as comparison with 
series expansions is possible. In both these limits we 
consider the "pure" t — J model which corresponds to 
t' = t" = and also a case with nonzero t' , t" . 



III. BOND-OPERATOR MEAN-FIELD THEORY 

The magnetic excitations in the magnetically disor- 
dered phase are triplons. To describe the triplons we em- 



ploy the spin-bond operator mean field technique. This 
approach has been previously applied to quantum dis- 
ordered systems such as bilayer antiferromagnets, spin 
chains, spin ladders, Kondo insulators etc ^ 44 i 45 i 53 ~ — . It 
is know n 44 ! 45 that this simple technique gives a QCP 
at (J±/J) c ~ 2.31, which is close to the exact value 
(J±/J) c — 2.525 known from Quantum Monte Carlo 
calculation s 46 ! 47 , series expansions 5 -^, and more involved 
analytical calculations 57 . 

In the present work we consider only the magnetically 
disordered phase since it proved to be sufficient to locate 
the LP. In principle one could extend the present study 
to the magnetically ordered phase. The bond-operator 
mean-field technique has been considered for the ordered 
phase of the bilayer antiferromagnet by both Normand 5 -^ 
and Vojta 3 . This could be considered in a future study. 
However, this is beyond the scope of the present work. 

All necessary equations describing the triplon dynam- 
ics have been derived previously^ 5 -. We present here a 
brief overview of the technique in application to the bi- 
layer antiferromagnet. Our aim is to establish a reliable 
description and to determine parameters which will be 
subsequently used in the SCBA calculations of a dressed 
hole, as addressed in Sec. [Vj One can certainly employ 
a more accurate Brueckner technique 5 - 7 -. However, this 
technique is more involved while the bond operator mean 
field approach has sufficient accuracy for our purposes 
and we chose it for simplicity. 

The bond-operator representation describes the system 
in a base of pairs of coupled spins on a rung, which can 
either be in a singlet or triplet (triplon) state: 

I s ) t = 4 1 o), = A=(\nh- nth) 
i t x )i = 4 1 o), = 4(itt)i- \uh) 

I t y )i = t\ y | >,= -^(in);+ \U)i) 

i t x )i = i i o >i= ^=(\n)i+ \m (2) 

where the four types of bosons obey the bosonic commu- 
tation relations. To restrict the physical states to either 
singlet or triplet, the above operators are subjected to 
the constraint, 

In terms of these bosons, the spin operators in each layer 
and can be expressed as 

s£ 2) = \{±*\ti« ± 4* - i E wl^) (4) 

where a represents respectively the components along the 
x, y, and z axes. Substituting the bond-operator repre- 
sentation of spins defined in Eq. dU into the Hj,j> in Eq. 
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(H|) we obtain 
Hj,j± = Hi + H2 



Hi 

#2 



H3 



7 Si Si 



H 4 
1 



Here the lattice spacing has been taken to be unity. The 
Bogoliubov coefficients u q and v q are given by 



t\ a tia) 



2n a 



i y^(s\s^ti a tj a + s i s jtiat] a + H.C.) 



(ii) 



The parameters /i and 5 are determined by the saddle 
point equations: 



(ho) 



(5) 



OH 



MFA 







dH 



MFA 



ds 







(12) 



The main issue is how to account for the hard core 
constraint (j3j). In principle this can be done via an infi- 
nite on-site repulsion of triplon excitations, however this 
technique is quite involved^. Our aim is the hole dis- 
persion, therefore for magnetic excitations the simplest 
possible technique which reproduces the magnetic QCP 
is sufficient. This is why we employ the simple mean 
field approac h 44 ! 45 that accounts the constraint (|3|) via a 
Lagrange multiplier \i in the Hamiltonian 

Hi -± Y,(4^ + tlji* ~ 1) 



Following Ref. |44], it is convenient to introduce the di- 
mensionless parameter d 



d- 



2Js 2 



(13) 



which results in the following self-consistency equations 



J_ 



N 



£ 



(6) 



^ ^ 1 ^ 

Jj_ £(-7 s J s i + jtiJia) ~ m£(>1 S ; + ^ioti* ~ l ) ■ 



V 1 + 

5 3 ^ 1 + dj q 
2 ~ 27V^T7T 



V 1 + 2d ^ 



3J± 



3J 

TV 



£ 



lq 



Further analysis is straightforward. We neglect the cu- 
bic and quartic Hamiltonians H3 and H4] replace singlet 
operators by numbers, (sj) = (si) = s (Bose-Einstein 
condensation of spin singlets); and finally diagonalize 
the quadratic in t Hamiltonian by performing the usual 
Fourier and Bogoliubov transformations 



V 1 + 2d ^ 



(14) 



Once the parameter d is determined from the first equa- 
tion, the values of s and \i can be calculated using the 
second/third relations. These parameters Eq. ([T4]) are 
then used to determine the excitation spectrum of the 
system 



£< 



£ 



Pa 
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a q L qa 



e ik ' ri a ko 



Vqt — q a 



(7) 

(8) 



J± 
4 



1 + 2d^ q 



1/2 



(15) 



Here N is the number of dimers; the diagonalized Hamil- 
tonian reads 



As a characteristic parameter of the magnetically dis- 
ordered ground state, we will be most interested in the 
value of the spin gap, i.e. the minimum energy gap of 
the triplon excitations, which is given by 



H M fa(v,s) = N 



3J ± s 2 



lis 



Ie^ 



J± 



1 - 2d 



1/2 



(16) 



where 



n q = 




A q = 


J± 
4 




Js 2 


B q = 





(cos(q x ) + cos(q y )) 



(10) 



(9) The results of numerical calculations of the spin gap is 
presented in Fig. [2] As J± decreases from infinity, the 
triplet gap at the wave vector q = (7r,7r) decreases, and 
the gap vanishes linearly at J±/J= (J_l/J) c , signalling 
the magnetic instability of the dimer phase. The critical 
value can be obtained by setting d = 1/2 in Eq. (p3J), 
which yields (Jj_/J) c ~ 2.3, in reasonable agreement 
with accurate Quantum Monte Carlo and series expan- 
sions result, (Jj_/J) c ~ 2.542a£L5£ Typical values of the 
parameters are presented in Table [III 
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FIG. 2: (color online) Triplon gap A at (tt,tt) as a function 
of J±/J. The solid and dashed areas along the J^/J-axis 
indicate the different topologies of the hole dispersion for the 
parameters t = 3.1, t' = —0.8, t" = 0.7, same as in Fig. 
EJb)-(d). Here the "4 hole pockets" region describes the four 
disconnected hole pockets with minima at k = (±7r/2, ±7r/2) 
while the " 1 hole pocket" region describes a single hole pocket 
with minimum at k = (tt, tt). 



J±/J 


s 


fi/J 




/ J A (7r;7r) /J 


2.31 


0.906 


-2.703 


-0.971 


0.021 


2.50 


0.918 


-2.756 


-0.881 


0.192 


2.70 


0.929 


-2.821 


-0.796 


0.384 


3.00 


0.943 


-2.938 


-0.688 


0.689 


3.50 


0.960 


-3.178 


-0.553 


1.218 


4.00 


0.971 


-3.459 


-0.459 


1.754 



TABLE II: Self consistent solution of the mean field param- 
eters as a function of J± at zero temperature. Calculations 
were performed by self consistently solving Eq. ([14)) and sub- 
stituting the resulting parameters in Eq. (|16[) . 



As one should expect the simple mean field approach 
described in this section does not give the correct critical 
index for the gap, see Fig. [2j however, for calculation of 
the hole Green's function the obtained accuracy is suf- 
ficient. One can think that the value of the "chemical 
potential" \i is a measure of correlations. Values of \i 
presented in Table [II] are surprisingly large. However, 
looking at Eq.© we see that | J± is just a redefinition 
of the ground state energy. Therefore, the true measure 
of correlations is fi + | J± . This quantity also presented 
in Table [II] is not too large and approaches zero when 
Jj_ — > oo. 



IV. HOLE-TRIPLON VERTEX 

At half filling (one electron per site) the model under 
consideration is equivalent to a Heisenberg model; specif- 
ically, it represents a Mott insulator with no long range 



antiferromagnetic order since we are considering the mag- 
netically disordered phase. In the triplon technique used 
in Section III the magnetically disordered state is built 
on the perturbative ground state consisting of isolated 
spin singlets (Bose condensation of spin singlets). Triplon 
virtual excitations accounted via Bogoliubov's transfor- 
mation gives the physical ground state in terms of the 
perturbative ground state. Since we will be using a di- 
agrammatic technique we need to define a hole in terms 
of the same perturbative ground state. Hence we define 
the hole creation operator with spin a by its action on 
the spin singlet bond i by a\ a 



«! t I s) =i o t) 

a\ l \s)=\oi) 



0) 
o), 



(17) 



where | 0) is vacuum. Note that the hole is created only 
at the upper plane as we have discussed in Section II. The 
electron annihilation operator in the upper plane can be 
expressed in terms of aj^, see e.g. Ref, 53 



i 

V2 



(p a Si + f; x ) + a\ a (pfftix + itiy) 



, (18) 



where p a = +(—), a =| (t) for a =j- (|). We also have 
to modify the hard core constraint (j3]) 



4 s i + Y ^iotioL + Y a L a i* = 1 



(19) 



however, this modification is not important for the single 
hole problem as well as it is not important for the low 
density hole problem. Note that the operators a\ a are 
required to obey fermionic anticommutation relations. 
This is important even in the single hole problem, since 
due to statistics, the sign of the hole dispersion is oppo- 
site to that of an electron. 

Substitution of (fT8|) in H t ,t',t" denned in (pQ) gives the 
Hamiltonian expressed in terms of hole and triplon oper- 
ators 

t- -s 2 

H t ,t',t»= Y K^> +a L a *^) 

JS ^ n . .+ r . _^ . JS , 

T 



2 ^ 



Js - 



Yl ^ijiTjj • (t\ a X tj a ) + T iy j ■ (tj a X ti a )) 

(iJ) 

J - Y KTi,i ' (t] a >< + Tjj ■ (4 X t^)) • 



Here 



(20) 



(21) 



7 



where r is the Pauli matrix. Eq. ([20]) contains hopping 
without changing the spin background (direct hopping), 
as well as spin-fluctuation assisted hopping, and exchange 
processes where the hole remains on its dimer. Consid- 
ering the direct hopping terms in Eq. ([20]) one obtains 
the bare hole dispersion 

ej™ = 2ts 2 lk + 2t'sS k + Zf'sSl , (22) 

where 

Ik = ^(cos(k x ) + cos{k y )) 
7 fc = cos(k x )cos(k y ) 

7* = \(cos(2k x ) + cos(2k y )) (23) 

The most important interaction is the hole-triplon ver- 
tex described by the effective Hamiltonian 

E E ' (4 +g <^^) + # - c - (24) 

feq era 

The interaction is due to the third and fourth lines in 
Eq. ([20]) . Taking matrix elements of these lines we find 
the verticies 



9k,q 


— 9k,q + 9k,q + 


9k,q + 9k,q 




2ts , 


+ Ik+qVq) 


9k,q 


= "^(^ 


9k,q 


= l=\lk U Q 


+ 7fe+g^) 


9k,q 


2t"s „ 

= {lk U Q 




9k,q 


Js 


+ 0*) 



Note, that on approaching the QCP the verticies diverge 
at q — )• (7r,7r). The verticies ([25]) differ both in the coef- 
ficients and in the kinematic structure from the verticies 
obtained in Ref. |43|. We believe that their verticies^ are 
not correct. 

The Hamiltonian ([20]) generates hole-double-triplon 
verticies as well. They are determined by the second, 5th 
and 6th lines of ([20]) . However, these verticies are due to 
quantum fluctuations of triplons and therefore they are 
very small. We neglect such verticies. 

Thus the effective Hamiltonian used in the next section 
for the calculation of the hole Green's function is of the 
following form 

H = E ^ + E 4° ] *L"k° (26) 

q k,a 
k,q,a,a 



V. HOLE GREEN'S FUNCTION 

The retarded Green functions for a hole is defined in 
the standard way 

POO 

G a (k,uj) = -i / (gs | a ka (t)al a (0) \ gs)e iut dt (27) 
Jo 

where | gs) is the ground state of the system. To describe 
the hole dressed by triplons we use the self-consistent 
Born approximation (SCBA) which disregards vertex 
corrections in Dysons' equation for the hole Green's func- 
tion; it amounts to an infinite order summation of non- 
crossing triplon loops (rainbow diagrams) in the hole self- 
energy, as presented diagramatically in Fig. [3] Vertex 
corrections shown in Fig. |4] are not accounted for in the 
SCBA. This approximation gives the leading term in the 
1/N expansion for the O(N) group^ 



+ 




FIG. 3: Dysons' equation in the the self- consistent Born ap- 
proximation. 




FIG. 4: Vertex corrections not included in the SCBA. 

This approximation has been widely used to study hole 
dynamics in the AF background^— . In the case of the 
AF background the vertex correction, Fig. HJ is equal to 
zero due to kinematic constraint o 18 i 19 i 6Q and hence the 
SCBA is very accurate. In the present case of the mag- 
netically disordered background the single loop vertex 
correction is nonzero. However, the correction is sup- 
pressed by the parameter 1/iV, where N = 3 is the num- 
ber of the triplon components^. To confirm the accuracy 
of the SCBA we compare results with that of numeri- 
cally exact dimer series expansions. Note that here we 
are working in terms of the true spin of the hole, while in 
the case of the AF background one has to work in terms 
of pseudospin. 

The hole dispersion in the magnetically ordered phase 
is strongly constrained by symmetry stemming from the 
magnetic ordering^!. However, here we consider the mag- 
netically disordered phase and hence the dispersion is free 
of such a constraint. The hole dispersion is purely deter- 
mined by the bare hole hopping and the underlying spin 
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dynamics, similar to the cuprates at high doping levels 
where the magnetic ordering becomes fully dynamic. 
The zero approximation Green's function is 



G°(k,o,) 



1 



(28) 



where is bare dispersion ([22]) . Summation of dia- 
grams Fig. [3] leads to the following Dyson's equation 

G a (k,uj) = ( cj-e^ 0) -3^^k-q, q ^(k - q,u-Q q )+iri 

^ q 

(29) 

where the factor 3 comes from three different polariza- 
tions of the intermediate triplon. Similar to the argument 
in Ref. [2l|, at infinite uj the Green's function satisfies 
G a = G®, therefore the sum rule 



1 f°° 

—Im / Gcr(k, uo)duo = 1 

* J-oo 



(30) 



is recovered. 

We solve Dyson's equation, Eq. ([29]) , numerically on a 
128 x 128 cluster with energy resolution Auj = 0.02 J and 
an artificial broadening r] = 0.02 J. The spectral function 
of a dressed hole is calculated by 



A(k,w) 



1 



Im[G(k,uj)} 



(31) 



VI. DIMER SERIES EXPANSIONS 

To check the accuracy of the SCBA approach we have 
also used series expansion method a 62 i 63 to compute the 
hole excitation energies. The approach we have used is to 
treat the J± terms exactly ('dimer expansion') and the 
other terms as perturbations. A linked-cluster expansion 
is used, wherein an effective Hamiltonian is derived for 
the 1-hole sector of each cluster. This yields a set of 
transition amplitudes for hopping of the hole through 
different lattice vectors, which can be summed to yield 
series for excitation energies for any wavevector k. The 
series are then analysed by standard Pade approximant 
methods. We refer the reader to Ref. [63| for further 
details of the method. In the present calculation, due to 
the presence of the further neighbor hopping terms t',t" 
the number of clusters grows rather rapidly. We have 
carried out the calculations to order 7, involving a total 
of 44312 distinct clusters with up to 7 dimers. 

Despite the rather short series, convergence is good 
provided the ^-parameters are not too large. By its na- 
ture, the series method will be accurate at small t-values 
whereas the SCBA is expected to be most accurate for 
large t. 




FIG. 5: (color online) The hole spectral functions at t = 
0.5, t' = t" = 0, for J ± = 2.31 (a) and J± = 3 (b). Here 
we show k = (f ,f) (black, solid), (7r,7r) (red, dashed), and 
(0,0) (blue, dotted). There is a significant distant incoherent 



part, the quasiparticle residues for J± = 3 are Z { 



: o.; 



s (o,o) 



0.32. 



(f>i) 



0.66, 



VII. NUMERICAL RESULTS 

First of all we remind the well known result^ - — that 
in the magnetically ordered phase at J± = the hole 
dispersion minimum is located at k = (±7r/2, ±7r/2). 
On the other hand at J± ^> J the hole dispersion is 
given by Eq. ([22j) with s = 1. We always take t > t',t", 
therefore the minimum of the dispersion is located at 
k = (±7r, ±7r). This implies that there is at least one 
Lifshitz point (LP) between J± = and J± = oo. Below 
we present and discuss our results and locate the LP. 

We start from the case of small hopping, t = 0.5, 
t' = t" = 0. At this small hopping we have two numeri- 
cal methods, the SCBA and dimer series expansions. We 
can compare results and hence judge the accuracy of the 
methods. This is especially important since vertex cor- 
rections are neglected in the SCBA. We plot in Fig. [5]the 
spectral functions at values J± = 2.31 and J± = 3. Note 
that J± = 2.31 is exactly the position of the QCP ob- 
tained from the mean field triplon analysis. We already 
pointed out that this position of the QCP is somewhat 
lower than that known from exact numerical calculations, 
J± = 2.52J. Since in the SCBA we use triplon spectra 
from Section Hill we must refer to the consistent position 
of the QCP. 

The spectra presented in FigJSJa) do not show any 
quasiparticle peaks but instead only power cuts are ob- 
served. This behaviour of the mobile hole spectral func- 
tion is similar to the Green's function of an immobile 
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(0,0) (f,f) (n,n) 

FIG. 6: (color online) The quasiparticle dispersion 6k along 
the nodal direction. For a convenient comparison we choose 
€7r,ir = 0. Calculations are performed for t = 0.5J, t' = 0.0, 
t" = 0.1J and for various values of the inter layer coupling 
J± . Points show results of the dimer series expansion method 
and curves show results of SCBA. Here we take J± — 2.31J 
(bottom), 2. 5 J, 2. 7 J, 3 J and 4 J (top). 



magnetic impurity at the QC P 59 i 64 , implying that re- 
markably the hole mobility does not influence this be- 
haviour. The power cuts imply that the spin is dis- 
tributed around the hole in a diverging cloud indicating 
SCS at the QCP (this will be discussed in more detail 
below). On the other hand spectra in Fig. [5fb) show 
quasiparticle peaks separated by the triplon gap A from 
the incoherent spectra. Figure [5(b) shows the dispersion 
minimum at k = (7r,7r), while the cut position in Fig. 
[5(a) is essentially the same for all momenta. Hence, we 
conclude that the position of the LP for these parameters 
coincides with that of the QCP. 

We also performed dimer series expansion calculations 
for both t = 0.5, if = t" = and t = 0.5J, if = 0.0, 
t" = 0.1 and compared with the results of the SCBA cal- 
culations. Note that the series expansion method allows 
one to determine only the quasiparticle dispersion. Nat- 
urally the method does not converge close to the QCP, 
as there are no quasiparticles there. However, at J± = 3 
the method works well and agreement between the SCBA 
and series is good, for example the SCBA band width is 
e(o,o) — e (7r,7r) = 0-40 and the series band width is 0.41 for 
the parameter set t = 0.5, t' = t" = 0. 

In Fig. [6] we show the comparison of the SCBA (lines) 
and dimer series expansions (points) for t = 0.5 J, t' = 
0.0, t" = 0.1. The series are well converged, even at 
the smallest value of J± . Hence the very good agreement 
between the two methods seen in Fig. [6]is strong evidence 
for the reliability of the SCBA approach for the present 
problem. In Fig. [6] the dispersion minimum at the QCP 
is k = (7r/2,7r/2). At higher J± the minimum moves 
towards k = (tt, tt) reaching k = (tt, tt) at about J± w 
4.5J, this is the LP. This is confirmed by the spectral 
functions and FS maps shown in Fig. [71 Thus the LP 
in this case is well separated from the QCP deep in the 
magnetically disordered phase. 

In the strong coupling limit t = 3.1 we rely on the 




FIG. 7: (color online) The hole spectral function at t — 0.5, 
t' = 0.0, t" = 0.1 for J ± = 2.31(a), 3.00(b) and 4.00(c). 
Values of momentum are k = (f , f ) (black, solid), (tt, tt) (red, 
dashed), and (7r,0) (blue, dotted). The quasiparticle residues 
for J± = 3 are ^(f ,f ) 0.63, Z (7T)7r) 0.90, Z (0>0 ) 0.32. 
The insets in Fig.(a)-(c) show maps of the hole dispersion. 
The dark region is minimum of the dispersion. 



SCBA since the series expansion does not converge in 
this strong coupling limit. Spectral functions for t = 3.1, 
t' = t" = are shown in Fig. [8] for J ± = 2.31 and J± = 
3.00. In this case there is no LP in the disordered phase 
since the spectral functions show the bottom of the band 
is always at k = (tt^tt). Hence for these parameters the 
LP is inside the magnetically ordered phase in agreement 
with Ref.— . It is worth noting that at the band bottom 
there are well defined quasiparticles even at the QCP, 
indicating no SCS for this choice of parameters. The 
spectra at k = (|, ^) and k = (tt, 0) have cuts at the 
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FIG. 8: (color online) The hole spectral function at t = 3.1, 
t' = t" = 0, for J ± = 2.31 (a) and J± = 3 (6). Values of the 
momentum are k = (77, 77) (black, solid), (7r,7r) (red, dashed), 
and (7r,0) (blue, dotted). The quasiparticle residues for J± = 
3 are = 0.31, Z^^) = 0.80, ^(^,0) — 0.35. The inset in 

Fig. (a) shows spectral functions in the broader energy range 
and the inset in Fig.(b) shows the map of the hole dispersion. 
The dark region is minimum of the dispersion. 



QCP, however, these are the high energy states which 
are irrelevant at small doping. 

The last and the most important set of parameters, 
t = 3.1, t' = —0.8, t" = 0.7, roughly corresponds to the 
parameters of the cuprate o 50 i 52 . Although the current 
study does not intend to simulate the quantitative de- 
tails of the cuprates, we choose the parameters relevant 
to the cuprates such that the interplay of various mech- 
anisms can be compared on a similar energy scale. The 
spectral functions are shown in Fig. [9] for several values 
of Jj_. The dispersion maps shown in the insets clearly 
demonstrate that the LP is located at J± « 3 well within 
the magnetically disordered phase. In Fig. [10] we also 
present plots of the quasiparticle dispersion and quasi- 
particle residue, which also clearly demonstrate that the 
LP is located at J± ~ 3. This topological transition is 
caused by fully dynamic antiferromagnetic correlations in 
the absence of any static magnetic order. This demon- 
stration of the possibility of the fully dynamic scenario is 
the first major conclusion of the present work. 

It is important to ask how sufficiently small a dop- 
ing is necessary for our analysis to be correct. In the 
present work we assume small doping, although it is pos- 



2. 
1.5 

< 

0.5 



2. 
1.5 
1. 

< 

0.5 





(c) 




(d) j 


: 

L 




r 








71 




• 


mm 





, A , J 




n 7 


— 71 

r 


-8 -7 


-6 -5 


-4 -3 





FIG. 9: (color online) The hole spectral function at t — 3.1, 
t' = -0.8, t" = 0.7 for Ji = 2.31(a), 2.50(b), 3.00(c) and 
4.00(d). Values of momentum are k = (f , f ) (black, solid), 
(7r,7r) (red, dashed), and (7r,0) (blue, dotted). The quasipar- 



ticle residues for J± = 3 are Zcjl il\ — 0.29, Z { 



^(7T,7T) 



0.77, 



0.20; The inset in Fig. (a) shows spectral functions in 
the broader energy range and the insets in Fig.(b)-(d) show 
maps of the hole dispersion. The dark region is minimum of 
the dispersion. 
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FIG. 10: (color online) Quasiparticle dispersion e& and residue 
Zk for different values of the interlayer coupling J± . Calcula- 
tions are performed with t = 3.1J, t' = — 0.8J and t" = 0.7 J. 
Here we considered J± = 2.5 J (red, solid), 3 J (blue, dashed), 
3. 5 J (black, dotted) and 4 J (blue, dot-dashed). 



sible to make a more quantitative estimate. It can be 
seen that practically we need x < 0.1 when the FS built 
on maps Fig|9](b) and Fig. 0(c) are topologically differ- 
ent. Another important question is why have the longer 
range hoppings have qualitatively changed the situation? 
The reason is this: at t' = t" = the LP is in the or- 
dered phase but already close to the QCP. A hopping 
t" > pushes the bare dispersion ([22]) at the nodal point 
k = ( § , § ) down helping magnetic fluctuations to form 
a small pocket. A pretty small positive t" is sufficient 
to shift the LP to the disordered phase. The role of t' 
is less important. The shift of the LP is due to the tun- 
ing of the longer range hoppings. The "tuning" has been 
performed by nature in the cuprates where the qualita- 
tive importance of t" is well known. These parameters 
give asymmetry between the hole and the electron dop- 
ing. Holes go to the nodal points while electrons go to the 
antinodal ones resulting in dramatically different Fermi 
surfaces and magnetic properties. We follow nature and 
rely on the same mechanism. 



VIII. SPIN-CHARGE SEPARATION 

The second major conclusion of the present work is 
related to the first one and concerns SCS at the QCP. 
Looking at Figj9ja), which is at the QCP, we observe 
that the lowest energy spectral function corresponding to 
k = (§,§), does not have a pole, but only a cut. Accord- 
ing to previous studies for the case of an immobile impu- 
rity, the presence of a cut indicates that the spin density 
is distributed in a power cloud around the hole 59,64 . Be- 
cause the Green's function in the present case is similar 
to that of the immobile impurity we directly project the 
results of Refs^^ to understand SCS at the QCP in 
the present case. When approaching the QCP from the 
magnetically disordered phase the quasiparticle residue 
approaches zero, Z oc A 2 , z « 0.4, as the triplon gap A 
approaches zero. The fraction of spin localized at hole 
goes to zero ex Z . The rest of spin is distributed around 
the hole over a disk of radius R oc 1/A. At r <C i? the 
spin density is oc l/r a , 1 < a < 1.5. Therefore the av- 
erage radius of the spin cloud (r) ~ R diverges at the 
QCP, indicating SCS. As one moves away from the QCP 
into the magnetically disordered phase, the quasiparticle 
peak appears, but the triplon gap is still significant as 
long as the system is close to the QCP (as can be seen 
in Fig. Efb)). This indicates a significant amount of spin 
is still relatively delocalized from the hole, and a rather 
smooth crossover to the spin charge recombined region. 
On the other side of the QCP, deep inside the AF phase, 
the hole interaction with a magnetic field is described 
by pseudospi n 4Q i 41 and this interaction implies a partial 
SCS^£. Evolution of the spin cloud when approaching 
the QCP from the AF phase is not clear at present. 

Because of the diverging size of the magnon cloud the 
hole effective mass also diverges at the QCP. Drawing 
analogy with the cuprates we note that the effective mass 
measured in MQO 28 diverges on approaching the QCP 
as identified by neutron scattering studie o 33 i 38 i 39 . 



IX. CONCLUSIONS 

We provide a detailed theory for the results presented 
in Ref. |42| where we considered the bilayer t — t' — t" — J 
model with strong interlayer coupling J±. In our anal- 
ysis we assumed that the hole doping is low and does 
not influence the magnetic dynamics and hence fill rigid 
bands formed by quantum magnetic fluctuations. We 
have argued that our results should be valid for dop- 
ing x < 0.10. Therefore, the magnetic dynamics are 
driven only by the interlayer coupling J± with a mag- 
netic QCP at J± ~ 2. 5 J separating the magnetically 
ordered and magnetically disordered phases. Compari- 
son between the hole dynamics described by the SCBA 
and series expansion methods are made explicitly, which 
shows excellent agreement. 

At J± — >• oo the hole Fermi surface is connected 
and centered at k = (7r,7r). We found that at a cer- 
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tain J^ p the AF correlations reconstruct the connected 
Fermi surface into four seperate pockets centered close 
to k = (±^,±7p. We have demonstrated that the LP 
where the topology of the Fermi surface is changed can be 
located within the magnetically disordered phase. The 
LP is driven by purely dynamic antiferromagnetic cor- 
relations in absence of any static magnetic order. The 
precise position of the LP at which this topological tran- 
sition occurs depends on hole hopping integrals t and t . 
This is because the reason behind the shifting of the hole 
pockets is due to magnetic fluctuations that diminish the 
nearest neighbor hopping t. The t and t then have the 
tendency to create hole pockets near k = (7r/2, 7r, /2). 
The precise position of LP depends on the balance of 
these two mechanisms. In fact, for certain values of t 
and t one may not see such transition if they are too 
weak to overcome the residual t. Nevertheless, for the 
parameters that are closest to the real cuprates, one sees 
such a transition at J±/J ~ 3, as indicated in Fig. [2j 

We have also demonstrated that if the LP is located in 
the magnetically disordered phase then, on approaching 
the magnetic QCP, the hole spin and charge separate. 
The separation scale is equal to the magnetic correlation 



length which diverges at the QCP. 

We are not aware of any 2D materials which have a 
low concentration of charge carriers and a magnetic QCP 
driven by a separation parameter. Perhaps such mate- 
rials will be synthesised in the future and/or the model 
can be realised in cold atoms. However, the most impor- 
tant outcome of the analysis is the demonstration of the 
principal possibility of the dynamic magnetic fluctuation 
driven topological transition and spin-charge separation. 
We anticipate that this calculation, although it does not 
describe the cuprate physics directly, indicates several 
important ingredients in the underdoped cuprates that 
are very likely to be the origin of FS reconstruction and 
phenemena related to SCS. 



X. ACKNOWLEDGEMENTS 

Numerical calculations performed in this work were 
done at the National Computational Infrastructure at the 
ANU, under project u66. 



1 S. I. Tomonaga, Prog. Theor. Phys. 5, 544 (1950). 

2 J. M. Luttinger, J. Math. Phys. 4, 1154 (1963). 

3 M. Vojta and K. W. Becker, Phys. Rev. B 60, 15201 
(1999). 

4 W. O. Putikka, R. L. Glenister, R. R. P. Singh, and H. 
Tsunetsugu, Phys. Rev. Lett. 73, 170 (1994). 

5 Y. C. Chen, A. Moreo, F. Ortolani, E. Dagotto, and T. K. 
Lee, Phys. Rev. B 50, 655 (1994). 

6 G. B. Martins, R. Eder, and E. Dagotto, Phys. Rev. B 60, 
R3716 (1999). 

7 X. G. Wen and P. A. Lee, Phys. Rev. Lett. 76, 503(1996). 

8 D. S. Marshall, D. S. Dessau, A. G. Loeser, C. H. Park, 
A. Y. Matsuura, J. N. Eckstein, I Bozovic, P. Fournier, A. 
Kapitulnik, W. E. Spicer and Z. X. Shen, Phys. Rev. Lett. 
76, 4841 (1996). 

9 X. J. Zhou, T. Cuk, T. De vereaux, N. Nagaosa, Z.-X. Shen, 
|arXiv:cond-mat/0604284| 

10 T Yoshida, X J Zhou, D H Lu, Seiki Komiya, Yoichi Ando, 
H Eisaki, T Kakeshita, S Uchida, Z Hussain, Z-X Shen and 
A Fujimori , J. Phys.: Condens. Matter 19, 125209 (2007). 

11 M. A. Hossain, J. D. F. Mottershead, D. Fournier, A. Bost- 
wick, J. L. McChesney, E. Rotenberg, R. Liang, W. N. 
Hardy, G. A. Sawatzky, I. S. Elfimov, D. A. Bonn and A. 
Damascelli, Nat. Phys. 4, 527 (2008). 

12 D. Fournier, G. Levy, Y. Pennec, J. L. McChesney, A. 
Bostwick, E. Rotenberg, R. Liang, W. N. Hardy, D. A. 
Bonn, I. S. Elfimov and A. Damascelli, Nat. Phys. 6, 905 
(2010). 

13 R-H He, X J Zhou, M Hashimoto, T Yoshida, K Tanaka, 
S-K Mo, T Sasagawa, N Mannella, W Meevasana, H Yao, 
M Fujita, T Adachi, S Komiya, S Uchida, Y Ando, F Zhou, 
Z X Zhao, A Fujimori, Y Koike, K Yamada, Z Hussain and 
Z-X Shen, New J. Phys. 13, 013031 (2011). 

14 H.-B. Yang, J. D. Rameau, Z.-H. Pan, G. D. Gu, P. D. 



Johnson, H. Claus, D. G. Hinks, and T. E. Kidd, Phys. 
Rev. Lett. 107, 047003 (2011). 

15 M. R. Norman, H. Ding, M. Randeria, J. C. Campuzano, 
T. Yokoya, T. Takeuchi, T. Takahashi, T. Mochiku, K. 
Kadowaki, P. Guptasarma and D. G. Hinks, Nature 392, 
157 (1998). 

16 S. Schmitt-Rink, C. M. Varma, and A. E. Ruckenstein, 
Phys. Rev. Lett. 60, 2793 (1988). 

17 C. L. Kane, P. A. Lee, and N. Read, Phys. Rev. B 39, 6880 
(1989) 

18 G. Martinez and P. Horsch, Phys. Rev. B 44, 317 (1991). 

19 Z. Liu and E. Manousakis, Phys. Rev. B 45, 2425 (1992). 

20 A. Ramsak, P. Horsch, and P. Fulde, Phys. Rev. B 46, 
14305 (1992). 

21 O. P. Sushkov, G. A. Sawatzky, R. Eder, and H. Eskes, 
Phys. Rev. B 56, 11769 (1997). 

22 I. M. Lifshitz Sov. Phys. JETP 11 (1960) 1130 [Zh. Eksp. 
Teor. Fiz. 38 (1960) 1569]. 

23 K. Y. Yang, T. M. Rice, and F. C. Zhang, Phys. Rev. B 
73, 174501 (2006). 

24 T. D. Stanescu and G. Kotliar, Phys. Rev. B 74, 125110 
(2006). 

25 N. Doiron-Leyraud, C. Proust, D. LeBoeuf, J. Levallois, J. 
B. Bonnemaison, R. Liang, D. A. Bonn, W. N. Hardy, and 
L. Taillefer, Nature 447, 565 (2007). 

26 B. Vignolle, A. Carrington, R. A. Cooper, M. M. J. French, 
A. P. Mackenzie, C. Jaudet, D. Vignolles, C. Proust and 
N. E. Hussey, Nature, 455, 952 (2008). 

27 N. Harrison, Phys. Rev. Lett. 102, 206405 (2009) 

28 J. Singleton, C. de la Cruz, R. D. McDonald, S. Li, M. Al- 
tarawneh, P. Goddard, I. Franke, D. Rickel, C. H. Mielke, 
X. Yao, and P. Dai , Phys. Rev. Lett, 104 086403 (2010). 

29 M. Le Tacon, G. Ghiringhelli, J. Chaloupka, M. Moretti 
Sala, V. Hinkov, M. W. Haverkort, M. Minola, M. Bakr, 



13 



K. J. Zhou, S. Blanco-Canosa, C. Monney, Y. T. Song, 
G. L. Sun, C. T. Lin, G. M. De Luca, M. Salluzzo, G. 
Khaliullin, T. Schmitt, L. Braicovich, and B. Keimer, Nat. 
Phys. 7, 725 (2011). 

30 K. Yamada, C. H. Lee, K. Kurahashi, J. Wada, S. Waki- 
moto, S. Ueki, H. Kimura, Y. Endoh, S. Hosoya, G. Shi- 
rane, R. J. Birgeneau, M. Greven, M. A. Kastner, and Y. 
J. Kim, Phys. Rev. B 57, 6165 (1998). 

31 V. Hinkov, S. Pailhs, P. Bourges, Y. Sidis, A. Ivanov, 
A. Kulakov, C. T. Lin, D. P. Chen, C. Bernhard and B. 
Keimer, Nature 430, 650 (2004). 

32 V. Hinkov, P. Bourges, S. Pailhs, Y. Sidis, A. Ivanov, C. D. 
Frost, T. G. Perring, C. T. Lin, D. P. Chen and B. Keimer, 
Nat. Phys. 3, 780 (2007). 

33 V. Hinkov, D. Haug, B. Fauqu, P. Bourges, Y. Sidis, A. 
Ivanov, C. Bernhard, C. T. Lin and B. Keimer, Science 
319, 597 (2008). 

34 M. H. Julien, A. Campana, A. Rigamonti, P. Carretta, F. 
Borsa, P. Kuhns, A. P. Reyes, W. G. Moulton, M. Horvatic, 
C. Berthier, A. Vietkin and A. Revcolevschi, Phys. Rev. B 
63, 144508 (2001). 

35 Y. Kobayashi, T. Miyashita, M. Ambai, T. Fukamachi and 
M. Sato, J. Phys. Soc. Jpn. 70, 1133 (2001). 

36 F. Coneri, S. Sanna, K. Zheng, J. Lord, and R. De Renzi, 
Phys. Rev. B 81, 104507 (2010). 

37 A. I. Milstein and O. P. Sushkov, Phys. Rev. B 78, 014501 
(2008). 

38 C. Stock, W. J. L. Buyers, Z. Yamani, Z. Tun, R. J. Bir- 
geneau, R. Liang, D. Bonn, and W. N. Hardy , Phys. Rev. 
B 77, 104513 (2008). 

39 D. Haug, V. Hinkov, Y. Sidis, P. Bourges, N. B. Chris- 
tensen, A. Ivanov, T. Keller, C. T. Lin and B. Keimer, 
New J. Phys., 12, 105006 (2010). 

40 S. A. Brazovskii and LA. Lukyanchuk, Zh. Eksp. Teor.Fiz. 
96, 2088 (1989) [Sov. Phys. JETP 69, 1180 (1989)]. 

41 R. R. Ramazashvili, Phys. Rev. Lett., 101, 137202 (2008); 
Phys. Rev. B 79, 184432 (2009). 

42 M. Holt, J. Oitmaa, W. Chen, and O. P. Sushkov, Phys. 
Rev. Lett. 109, 037001 (2012). 

43 C. Briinger and F. F. Assaad, Phys. Rev. B 74, 205107 
(2006). 



44 Y. Matsushita, M. P. Gelfand, and C. Ishii, J. Phys. Soc. 
Jpn. 68, 247 (1999). 

45 D. K. Yu, Q. Gu, H. T. Wang, and J. L. Shen, Phys. Rev. 
B 59, 111 (1999). 

46 A. W. Sandvik and D. J. Scalapino, Phys. Rev. Lett. 72, 
2777 (1994). 

47 A. W. Sandvik, A. V. Chubukov, and S. Sachdev, Phys. 
Rev. B 51, 16483 (1995). 

48 J. M. Tranquada, G. Shirane, B. Keimer, S. Shamoto, and 
M. Sato, Phys. Rev. B 40, 4503 (1989). 

49 D. Reznik, P. Bourges, H. F. Fong, L. P. Regnault, J. 
Bossy, C. Vettier, D. L. Milius, LA. Aksay and B. Keimer, 
Phys. Rev. B 53, R14741 (1996). 

so y Tokura, S. Koshihara, T. Arima, H. Takagi, S. Ishibashi, 
T. Ido, and S. Uchida, Phys. Rev. B 41, 11657 (1990). 

51 M. Greven, R. J. Birgeneau, Y. Endoh, M. A. Kastner, B. 
Keimer, M. Matsuda, G. Shirane, and T. R. Thurston , 
Phys. Rev. Lett. 72, 1096 (1994). 

52 O.K. Andersen, A.I. Liechtenstein, O. Jepsen, and F. 
Paulsen, J. Phys. Chem. Solids 56, 1573 (1995). 

53 C. Jurecka and W. Brenig, Phys. Rev. B 63, 094409 (2001). 

54 R. Eder, Phys. Rev. B 57, 12832 (1998). 

55 Y. Saito, A. Koga, and N. Kawakami, J. Phys. Soc. Jpn 
72, 1208 (2003). 

56 ZhengWeihong, Phys. Rev. B 55, 12267 (1997). 

57 V. N. Kotov, O. Sushkov, ZhengWeihong , and J. Oitmaa, 
Phys. Rev. Lett. 80, 5790 (1998). 

58 B. Normand and T. M. Rice, Phys. Rev. B 56, 8760 (1997). 

59 O. P. Sushkov, Phys. Rev. B 62, 12135 (2000). 

60 O. P. Sushkov, Phys. Rev. B 49, 1250 (1994). 

61 W. Chen, O. P. Sushkov and T. Tohyama, Phys. Rev. B 
84, 195125 (2011). 

62 J. Oitmaa, C. J. Hamer, and ZhengWeihong, Phys. Rev. 
B. 60, 16364 (1999) 

63 J. Oitmaa, C. J. Hamer and W. Zheng, Series Expansions 
for Strongly Interacting Lattice Models, (Cambridge Uni- 
versity Press, New York, 2006) 

64 M. Vojta, C. Buragohain, and S. Sachdev, Phys. Rev. B 
61, 15152 (2000). 



